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OPTIMAL  CONTROL  OF  UNSTEADY  STOKES  FLOW  AROUND  A  CYLINDER  AND 
THE  SENSOR/ACTUATOR  PLACEMENT  PROBLEM 

JOSIP  LONCARIC* 

Abstract 

Abstract.  Effective  placement  of  sensors  and  actuators  is  of  crucial  importance  in  flow  control.  Instead  of 
using  combinatorial  search  to  identify  optimal  locations,  we  pose  a  related  problem  of  polynomial  complexity. 

If  one  could  sense  everything  and  actuate  everywhere,  what  should  one  do?  Using  the  unsteady  2D  Stokes 
flow  around  a  cylinder  as  an  example,  we  obtain  the  analytic  solution  of  an  optimal  distributed  control 
problem  and  describe  its  spatial  structure.  At  low  circumferential  wavenumbers  or  close  to  the  cylinder  wall, 
boundary  vortex  generators  are  shown  to  be  more  effective  than  colocated  vorticity  damping.  This  analytic 
solution  has  also  been  used  to  test  numerical  methods,  demonstrating  the  importance  of  using  discretization 
which  resolves  all  eigenfunctions  of  interest. 

Key  words,  optimal  flow  control,  exterior  Stokes  flow,  sensor/actuator  placement,  design 

Subject  classification.  Applied  and  Numerical  Mathematics,  Controls 

1.  Introduction.  Given  the  task  of  controlling  a  physical  system  described  by  partial  differential  equa¬ 
tions,  one  must  first  select  sensors  and  actuators,  then  devise  a  controller  for  the  resulting  distributed  system. 
Although  intuition  and  experience  can  sometimes  help  in  the  design  stage,  this  approach  fails  when  prior 
experience  is  lacking.  Confronted  with  a  radically  different  problem,  such  as  the  possibility  of  using  several 
small  sensors  and  actuators,  one  quickly  finds  that  the  number  of  possible  configurations  grows  combinato- 
rially.  For  example,  [11]  investigated  the  design  optimization  problem  of  placing  M  =  8  active  struts  into  a 
structure  with  AT  =  1507  possible  locations,  resulting  in  a  discrete  search  space  of  approximately  6.5  X  10^0 
possible  designs.  The  search  space  grows  dramatically  as  either  A/'  or  M  is  increased. 

While  this  formulation  of  the  problem  may  be  natural  in  structures  composed  of  discrete  components 
such  as  struts,  a  different  approach  is  needed  in  continuum  problems  where  the  number  of  candidate  locations 
is  potentially  infinite.  We  shall  adopt  the  point  of  view  that  the  search  space  for  placement  of  sensors  and 
actuators  may  be  reduced  by  performing  a  thought  experiment  in  optimal  control.  This  thought  experiment 
asks  the  following  question:  If  one  could  sense  everything  and  control  everything,  what  should  one  do? 

While  distributed  control  is  typically  not  feasible,  it  leads  to  the  optimal  design  which  we  can  try  to 
approximate  by  feasible  designs.  We  anticipate  that  the  optimal  feedback  law  favors  some  spatial  regions 
over  others.  This  can  happen  if  the  system  dynamics  exhibit  spatially  localized  behavior,  such  as  boundary 
layers  in  flow  control  problems,  or  through  a  particular  choice  of  the  optimality  criterion.  In  either  case,  we 
espect  that  the  optimal  distributed  feedback  law  will  be  of  the  form 

K{a,  s)x{s)ds 

where  u{a)  is  the  optimal  control  input  at  each  actuation  location  a  in  O,  rr(5)  is  the  state  variable  at  sensing 
location  s,  «(a,  s)  is  the  optimal  feedback  kernel  and  the  integral  ranges  over  the  entire  volume  of  the  domain 
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fl.  As  observed  by  [1,  12],  this  integral  may  be  approximated  by  neglecting  the  pairs  of  locations  for  which 
K{a,  s)  is  small.  The  remaining  pairs  of  locations  where  «(a,  s)  is  large  define  a  reduced  design  search  space. 
The  thought  experiment  can  also  be  applied  recursively,  until  the  search  for  feasible  and  cost  effective  designs 
becomes  manageable. 

Prior  work  on  this  problem  focused  on  controllability  and  observability  for  distributed  parameter  sys¬ 
tems,  leading  to  to  the  concept  of  strategic  actuators  and  sensors  [5,  6].  This  property  can  be  completely 
discontinuous.  For  the  heat  equation  on  the  unit  interval,  a  point  actuator  or  sensor  is  strategic  if  and  only 
if  its  coordinate  is  an  irrational  number. 

The  computational  method  introduced  in  [10]  is  most  suitable  for  design  problems  where  a  small  set  of 
candidate  locations  is  already  heuristically  determined.  Combinatorial  growth  of  computational  complexity 
is  avoided  by  restricting  attention  to  candidate  sensor/actuator  pairs.  The  pairwise  dimensions  of  the 
intersection  of  controllable  and  observable  subspaces  are  computed  first,  and  pairs  are  ranked  in  terms  of 
effectiveness  and  versatility.  For  each  significant  vibration  mode  of  the  testbed  structure,  the  best  pairs  were 
chosen  based  on  a  combined  index  defined  in  terms  of  the  controllability  and  observability  Gramians. 

Computing  the  optimal  n{a,s)  in  (1.1)  also  leads  to  a  problem  of  polynomial  complexity,  but  without 
initial  assumptions  about  candidate  locations.  Given  a  discretization  of  space  into  locations,  the  linearized 
dynamics  can  be  described  by  a  iV^  x  matrix  and  the  optimal  feedback  found  in  at  most  0{N^)  steps. 
By  contrast,  in  placing  M  sensors  or  actuators  the  discrete  search  approach  would  require  the  evaluation  of 


(1.2) 
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possible  designs  when  N  =  100  and  M  =  10.  The  cost  of  computing  the  optimal  feedback  kernel  is  much 
lower,  of  order  =  10^^.  While  not  completely  impossible,  this  is  still  daunting  even  in  2D  problems, 
where  the  computational  cost  would  be  of  order  =  10^^.  Further  reduction  of  computational  complexity 
is  needed. 

Instead  of  reducing  spatial  resolution,  we  shall  focus  on  constructing  a  reduced  order  model  of  a  particular 
system.  We  shall  seek  to  construct  a  sequence  of  reduced  order  models  by  using  a  modal  approximation. 
This  approach  may  be  rigorously  justified  provided  that  the  system  belongs  to  the  class  of  spectral  systems 
defined  in  [2,  3].  In  fact,  we  shall  show  how  system  dynamics  of  the  unsteady  Stokes  flow  around  a  cylinder 
may  be  diagonalized,  derive  the  exact  optimal  feedback  kernel,  and  then  construct  its  analytic  approximation 
whose  worst  case  performance  loss  is  less  than  0.026  percent.  We  begin  by  formulating  the  problem. 


2.  Exterior  Stokes  flow  problem.  Much  of  our  understanding  of  the  full  Navier-Stokes  equations 
begins  with  the  Stokes  flow,  which  represents  their  linearization  around  the  motionless  state  [13].  Incom¬ 
pressible  fluid  dynamics  can  be  thought  of  as  the  response  of  the  Stokes  flow  to  forcing  by  the  nonlinear 
term.  While  the  nonlinear  term  conserves  the  kinetic  energy,  a  vorticity  perturbation  can  grow  via  the  vortex 
stretching  mechanism  in  3D.  In  2D,  this  mechanism  is  absent  and  the  vorticity  is  conserved  along  particle 
paths  when  viscosity  and  forcing  are  neglected.  Consequently,  the  total  enstrophy  |||a;|p  cannot  increase 
under  the  action  of  the  nonlinear  term  alone  [7].  Instead,  additional  vorticity  is  created  at  the  solid  bound¬ 
aries  in  a  process  governed  by  the  viscous  sublayer  and  the  Stokes  flow  equations.  A  correct  understanding 
of  the  boundary-fluid  interaction  becomes  particularly  significant  when  boundary  control  is  contemplated. 

The  flow  of  incompressible  viscous  fluid  around  a  cylinder  is  a  prototype  exterior  flow  problem.  For 
control  purposes  at  low  Reynolds  numbers,  vorticity  creation  at  the  solid  boundaries  and  viscous  dissipation 
are  the  two  physical  processes  of  primary  importance,  while  the  conservative  nonlinear  term  is  secondary. 
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This  suggests  that  a  simpler  model  in  which  the  nonlinear  term  is  neglected  would  still  capture  the  essential 
aspects  of  this  flow  control  problem. 

Consider  the  corresponding  2D  Stokes  flow  in  streamfunction  representation  described  by  the  equations 


(2.1) 

V  =  — z  X  V'lp 

(2.2) 

a;  =  Vxv  =  —A'lp 

BllJ 

(2.3) 
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where  ^  is  the  streamfunction  whose  gradient  rotated  by  -'7r/2  is  the  velocity  field  v.  The  boundary 
conditions  are 


(2.4) 

V'|r=l  =  0 

(2.5) 

=0 

r=l 

(2.6) 

lim  =  0 

r— >oo 

As  is  customary  in  2D  flows,  the  velocity  v  will  be  thought  of  as  a  2-vector  but  ^  and  the  vorticity  uj  will 
be  considered  scalars.  Moreover,  we  shall  interpret  v  as  a  finite  energy  flow  perturbation. 

At  this  point  it  is  useful  to  note  that  the  scalar  second  order  elliptic  equation  =  —a;  must  satisfy  more 
than  two  boundary  conditions,  which  is  only  possible  if  the  vorticity  uj  satisfies  certain  integral  compatibility 
conditions.  As  pointed  out  by  [9],  this  compatibility  is  achieved  by  the  creation  of  vortex  sheets  at  the  solid 
boundary.  This  vorticity  creation  process  at  the  wall  couples  to  the  exterior  flow  through  viscosity,  diffusing 
the  vortex  sheet  outwards  in  an  arbitrarily  short  time.  While  the  implicitly  created  vortex  sheets  involve 
delta  functions  which  present  a  challenge  to  numerical  approximations,  we  intend  to  consider  a  family  of 
reduced  order  models  which  approximate  the  full  dynamics  by  neglecting  the  highly  stable  modes.  As  we 
shall  show  later,  the  retained  modes  can  be  expressed  analytically  in  terms  of  smooth  functions  which  can 
be  approximated  numerically  without  much  difficulty. 

We  are  interested  in  flows  for  which  the  initial  vorticity  is  contained  within  the  computational  domain 
of  radius  R  and  the  diffusion  of  vorticity  across  the  artificial  boundary  v  =  Ris  negligible.  To  that  end,  let 
us  introduce  coordinate  variables  p  =  log(r)  and  then  apply  the  Fourier  transform  in  the  6  direction.  The 
relevant  operators  may  be  written  as  follows: 


This  factorization  of  the  Laplacian  in  Fourier  representation  separates  the  modes  which  grow  and  decay  as 
r  — ^  00.  In  complex  analysis,  this  is  equivalent  to  the  Wiener-Hopf  factorization,  which  has  been  used  to 
study  boundary  layer  growth  near  the  leading  edge  of  a  flat  plate  [8],  but  it  will  also  prove  useful  in  our  2D 
problem.  We  note  that  this  approach  generalizes  to  higher  dimensions  via  the  Calderon  projector  [15] . 

As  shown  in  [14]  for  the  problem  Au  =  /,  whenever  the  forcing  /  vanishes  outside  the  computational 
domain,  the  boundary  condition  u  — >  0  at  infinity  is  exactly  represented  by  the  countably  many  conditions 

(2.9)  ( \k\^  Ukip)  =0,  A:  =  0,±1,±2,... 

/  P=iog(ie) 

(2.10)  uo(log(it:))  =  0 
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applied  at  the  artificial  boundary  of  the  computational  domain.  These  boundary  conditions  are  nonlocal 
and  are  best  applied  in  the  Fourier  representation. 

3.  Pseudodiflferential  representation.  Let  us  also  introduce  the  following  pseudodifferential  repre¬ 
sentation. 

Definition  3.1.  The  pseudodifferential  representation  ^  of  the  flow  state  with  the  streamfunction  ^  is 
given  by 

(3.1)  6(p)  =  ■  +  l^l) 

This  definition  implies  that 

(3.2)  (p)  =  -vl{p)-i  sgn{k)vl  {p) 

where  and  are  the  circumferential  and  radial  components  of  the  solenoidal  velocity  field,  respectively. 
We  note  that  ^  is  a  real  scalar  field.  Moreover,  unlike  the  streamfunction  or  vorticity,  ^  scales  with  velocity 
as  the  physical  length  scale  is  changed. 

Theorem  3.2.  Given  a  function  ^kip)^  I'he  streamfunction  'ipkip) 


(3-3) 

Mp)=  re-l'''^"-")efe(<^)e"rf<7 
Jo 

(3.4) 

The  proof  follows  by  application  of  the  variation  of  constants  formula  and  the  fact  that  the  solid  boundary 
r  =  1  is  the  streamline  ^  =  0, 

Furthermore,  Vk{p)  —  0  clearly  implies  that  ^k{p)  ~  0.  The  converse  need  not  be  true  since  ^k{p)  =  0 
implies  only  that  u^(p)  =  — sgn(fc)i;^(p)  which  has  nonzero  solenoidal  solutions  derived  from  corner  flows 
=  1/rl^l  for  r  >  R.  However,  as  an  easy  consequence  of  the  above  theorem,  we  conclude  that  when  ^ 
vanishes  outside  a  bounded  domain,  v  decays  to  zero  at  infinity. 

Corollary  3.3.  Whenever  the  support  of  ^  is  contained  within  a  bounded  domain  r  <  R,  the  corre¬ 
sponding  flow  field  satisfies  the  boundary  condition  v  ^  0  at  infinity. 

We  can  now  prove  a  stronger  result  which  holds  even  when  ^  does  not  have  compact  support. 
Theorem  3.4.  The  following  boundary  conditions  on  ^  and  v  at  infinity  are  equivalent  for  continuous 
functions 

(3.5)  lim  £  =  0  lim  v  =  0. 

p—*oo  ^  p-*oo 

The  implication  <=  is  a  trivial  consequence  of  the  definition  of  £.  To  obtain  the  implication  =>,  note 
that  the  velocity  field  is  linear  in  £.  Let  us  write  £  =  £o  +  ^oo  where  the  £oo  vanishes  for  p  <  log{R)  and  £o 
vanishes  for  p  >  \og{R) .  The  corresponding  solenoidal  velocity  field  is  v  =  vq  +  Vo©  where  Vq  decays  to  zero 
at  infinity  by  the  previous  corollary.  For  Vq©  we  use  the  boundary  condition  -0  =  0  at  r  =  1  to  obtain  the 
Fourier  coefficients  of  the  radial  velocity 

(3.6)  vlo,kip)  = 

y  r 
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(3.7) 

(3.8) 

for  p  >  \og{R).  Since  ^  decays  to 


=  4  re~''''^'’-")e"eoo,fe(<T)dcT 

Jo 

=  ik  f 

J\og(R) 

at  infinity,  define  e{k^R)  by 


(3.9) 


e{k,R):=  sup  \^k{p)\ 

p>\og{R) 


which  leads  to  the  inequality 


(3.10) 


C,fc(p)l  <  e{k,R)\k\  [ 
J\c 


p 

log(H) 


Taking  the  limit  p  oo,  we  obtain 

(3.11)  Ji^  b^.fc(p)l  <  €(fc, -R) 


which  holds  for  all  fc,  R.  Since  ^  is  continuous,  its  Fourier  series  converges  absolutely  and  uniformly,  and  so 
does  the  Fourier  series  for  Given  that  ^  >  0  at  infinity,  the  conclusion  >  0  follows  by  taking  the 

limit  R^  oo.  Prom  the  definition  of  we  then  obtain  that  the  tangential  component  0  and  the  proof 
is  complete. 

At  this  point  the  main  advantage  of  the  pseudodifferential  representation  has  become  clear.  We  shall 
represent  flows  in  terms  of  where  ^  satisfies  Dirichlet  boundary  conditions.  By  contrast,  the  velocity- 
vorticity  and  the  streamfunction-vorticity  representations  involve  the  difficulty  of  determining  the  vorticity 
created  at  the  solid  boundary. 

While  the  pseudodifferential  representation  has  been  defined  in  terms  of  its  Fourier  coefficients,  an  equiv¬ 
alent  physical  interpretation  in  terms  of  singular  integral  operators  can  be  found.  The  resulting  equations 
explicitly  demonstrate  the  non-local  nature  of  the  pseudodifferential  representation. 

Theorem  3.5.  Given  a  smooth  streamfunction  the  pseudodifferential  representation  ^  is  given  in 
terms  of  physical  coordinates  by 


(3.12) 


r 


(  dipjp,  0) 

V  dp 


d'^-ipjp,  e-4)) 

de‘^ 


log 


The  proof  follows  by  applying  the  convolution  theorem  to  the  terms  in  the  integral  and  integrating  by 
parts,  since  the  function  —2  log  |  sin(a:/2)|  has  Fourier  coefficients  l/\k\  for  nonzero  k.  For  nonsmooth  -0,  this 
equality  may  still  be  interpreted  in  the  sense  of  distributions,  or  taken  as  an  alternative  definition  of  This 
result  can  also  be  expressed  in  terms  of  the  velocity  fields. 

Corollary  3.6.  The  pseudodifferential  representation  ^  is  related  to  smooth  divergence-free  velocity 
fields  by 


(3.13) 


1 

^ip,e)  =  -v%p,e)  +  -  / 
^  Jo 


dv''{p,  0  -  (t>) 


log 


.  0 

sin  — 
2 


d<f). 
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4.  Reformulated  Stokes  flow  problem.  Besides  the  field  we  define  the  following  three  operators 

(4.1)  Tlfip)  =  rf{p) 

(4.2)  Bik\fkip)  =  “  l^l)  Mp') 

(4.3)  ^\k\fkip)  =  fk{p) 

and  their  particular  inverses 

(4.4)  n-^fk{p)  =  ^ 

poo 

(4.5)  B^,\Mp)  =  -  /  e-l''l(-P)/,(<T)da 

Jp 

(4.6)  ^|fc|Vfc(p)  = 

Given  these  definitions,  the  factorization  (2.8)  of  the  Laplacian  may  be  compactly  written  as 

(4.7)  A|fc|/fc(p)  =  n-%k\T\k\fk{p). 


We  shall  also  be  using  the  commutativity  properties 

(4.8)  Biki+aT^""  = 

(4.9) 

(4.10)  BaTb  =  ThBa  . 


While  the  above  definitions  are  given  in  terms  of  Fourier  components,  physical  interpretation  of  these 
operators  can  be  found.  The  following  theorem  is  an  easy  consequence  of  the  convolution  theorem  once  we 
observe  that  for  positive  a  the  Fourier  coefficients  represent  the  function 


(4.11) 


def  sinh(cr) 
cosh(c7)  —  cos{(j)) 


for  (7  >  0 


whose  limit  as  a  ^  O"*"  has  all  Fourier  coefficients  equal  to  1  provided  that  it  is  defined  to  be 


(4.12) 


g{0,(l))  =  7rJ(sin((/>/2)). 


The  function  g{a,  0)  may  be  recognized  as  the  Poisson  kernel. 

Theorem  4.1.  The  above  operators  and  their  particular  inverses  have  the  following  physical  represen- 


tations: 

(4.13) 

TZ  :  fip,e)  >-^rf{p,e) 

(4.14) 

(4.15) 

^  m  9f{p,e)  1  p^d^fiP, 9 

dp  2,1  SS^ 

poo  1  p2'K 

(4.16) 

:  fip,  ~  f  f  ^^f^P  +  <^^9-‘f>)  d(pda 

(4.17) 

.r  .N  9f{p,e)  ,  1  r^d^f{p,9-(f>)^,^,,^ 

dp  +2,1 

pp  -1  p2'K 

•  f{p,9)>-^  ^  9i(^,‘f>)f{p-(^,9-^)d(t>da 

(4.18) 
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where  r  =  e^,  h{(j))  =  log  (4sin((^/2)^)  and  ^(cr,  </>)  is  defined  as  above. 

After  shortening  the  notation  by  making  the  dependence  on  p  implicit,  the  model  problem  is  described 
by  the  following  theorem: 

Theorem  4.2.  The  pseudodifferential  representation  of  the  Stokes  flow  around  a  cylinder  is  described 
by  the  following  equations: 


(4.19)  =  ^|fc| 

(4.20)  Wk  = 

(4.21)  vl  =  (^Ife|  -  e,,,)  V-fe  =  i 

(4.22)  vt  =  +  J-|,|)  V’fc  =  (lfc|^|^|+i  -  & 

(4.23)  ^  =  A|n.|_i|^fc 

with  Dirichlet  boundary  conditions  on 

Equation  (4.19)  follows  from  the  definition  of  ^  and  the  boundary  condition  ^  =  0  at  r  =  1.  Equations 
(4.21)  and  (4.22)  also  follow  immediately.  The  equation  (4,20)  is  verified  as  follows: 


(4.24) 

(4.25) 

(4.26) 

(4.27) 


Wfc  =  -Aipk 

—  -T^~^B\k\-l'^~^^\k\'>Pk 


This  result  shows  that  whenever  ^  vanishes  for  r  >  R,  so  does  u).  Conversely,  given  the  boundary  condition 
V  — »  0  at  infinity,  one  may  invert  the  equation  (4.20)  and  show  that  whenever  w  vanishes  for  r  >  R^,,  so 
does  £.  Both  properties  are  desirable  given  the  bounded  computational  domain. 

Assuming  that  the  initial  vorticity  distribution  is  contained  within  a  smaller  domain  r  <  R^  <  R, 
diffusion  of  vorticity  out  of  the  computational  domain  is  insignificant  for  t  <C  {R  —  Ru,)"^.  By  choosing 
R'>  Ru  +  VT  this  reformulated  problem  can  be  made  applicable  for  all  t  <  T.  Provided  that  the  initial 
disturbance  decays  to  almost  zero  by  the  time  t  =  T,  the  forcing  at  the  artificial  boundary  r  —  R  will  again 
be  negligible.  In  the  limit  R—yoo,  the  pseudodifferential  representation  satisfies  the  the  equivalent  boundary 
conditions  over  the  entire  exterior  domain. 

Finally,  equation  (4.23)  is  obtained  from  (2.3)  as  follows: 


(4.28) 

(4.29) 

(4.30) 

(4.31) 

(4.32) 


du> 

'dt 

d^k 


=  Auj 


Bi 


lfcl-l 


dt 

d^ 

dt 


dt 

djk 

dt 


where  the  last  step  follows  from  the  enforced  artificial  boundary  condition  ^  =  0  for  r  >  ii.  As  R  -y  oo, 
another  justification  of  this  step  follows  from  the  boundary  conditions  on  ^  and  the  fact  that  the  diffusion  of 


7 


any  bounded  compactly  supported  initial  vorticity  distribution  is  so  slow  that  T\k\uJk{p)  must  remain  zero 
at  infinity  over  all  finite  subsequent  time  intervals. 

We  can  also  interpret  the  relationship  between  -0,  ^  and  uj  in  terms  of  Cartesian  coordinates. 
Corollary  4.3.  Denoting  the  exterior  points  by  Cartesian  vectors  jc  ory,  the  following  identities  hold: 

■Jy“ 


(4.33) 

^(x)  =  ^  ^(y)  1 

2i7r  J  Jl<|y|<|x|  1, 

(4.34) 

f  1  ^(y)  j: 

(4.35) 

dip{x)  1  /  < 

(4.36) 

‘ 

^  ^  dr  r  27t 

y|‘ 


ds‘  2|x||yl 

f  <fay) 


log 


/|x||y|  -x-y 
I  2|x||y| 


ds 


where  r  =  |x|  and  along  the  contour  of  integration  y  =  y(s)  where  s  =  r9  is  the  distance. 

The  proof  is  based  on  theorems  4,1  and  4.2,  and  the  observation  that  in  Cartesian  coordinates  the 
functions  g  and  h  simplify  to 

ix|2  -  |xf 


(4.37) 

(4.38) 


g{p  -p',e-  e')  = 


|x 


—  v'|2 


h{e  -  0')  =  log 


|x||x'|  -x-x'\ 
2|x||x'|  ) 


along  |x|  =  |x'|. 


We  shall  consider  the  pseudoenergy  as  equivalent  to  the  kinetic  energy  |||v||^.  The  following 

theorem  shows  that  this  is  justified  whenever  V^(p,  ^)  approaches  a  function  of  p  alone  as  p  ^  6o.  This 
situation  applies  to  flows  of  interest,  whose  initial  vorticity  distribution  is  compactly  supported  and  total 
kinetic  energy  is  finite. 

Theorem  4.4.  The  total  pseudoenergy  is  equal  to  the  total  kinetic  energy  and  ||^|p  =  ||v|p  whenever 
the  following  sufficient  condition  holds: 


(4.39) 


p2Tr 

lim  / 


dip 


de  ^0 


Proof: 


(4.40) 

nOC 

lien"  =  / 

Jo 

>  n2TV 

/  dO  dp 

Jo 

(4.41) 

>  r2Tr 

Jo  ^ 

d^rp? 

dO  dp 

(4.42) 

=f 

'?( 

dipk 

dp 

^  +  k'^\i’k?  +  \k\(^ 

(4.43) 

-r 

)  ^27r 

i  ' 

dedp+  yjkWipkl'^ 

k 

(4.44) 

=  l|v|| 

2 

#fc  ,  7  dipk 


p=0 


dp 


since  at  each  p  the  infinite  sum  is  bounded  by 
(4.45) 


1  /‘Stt 

0  <  J2^k\\ipk\'^  <  =  ^ 

k  k 


dip 


de 


dd 
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which  vanishes  at  p  =  0  due  to  the  boundary  conditions  and  as  p  oo  due  the  assumption  of  the  theorem. 
Although  an  even  milder  assumption  would  suffice,  we  shall  be  satisfied  with  this  result  for  the  sake  of 
simplicity. 

5.  Diagonalization  of  the  ^  dynamics.  At  this  point,  we  can  solve  the  eigenvalue  problem  and 
obtain  the  modes  of  the  pseudodifferential  representation  of  the  Stokes  flow.  First  we  shall  consider  a 
bounded  computational  domain,  and  then  let  i?  ^  oo  to  include  the  full  exterior  domain. 

On  a  bounded  computational  domain,  system  dynamics  are  governed  by  the  operator  defined  by 

(5.1)  Akfk  = 

over  the  range  1  <r  <  R  and  0  otherwise. 


Z  k=0  m=0. 39, 0.73, 1.1  w  k=0  m=0 . 39 . 0 . 73 , 1 . 1 


Fig.  5.1.  Thjpical  modes  of^  (left)  andco  (right)  vs.  radius. 

Theorem  5,1.  The  normalized  solutions  of  the  equation  {Ak  4-  =  0  restricted  to  the  range 
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I  <r  <  R  are  linear  combinations  of  the  Bessel  functions 


(5.2) 


^nm(^)  — 


Yn{mr)Jn{m)  -  Jn{mr)Yn{m) 

+  y„(m)2 


where  n  |1A:|  -  1|  and  m  is  a  root  of  the  equation  ZnmO-)  =  Znm{R)  =  0- 

This  theorem  is  proved  by  direct  substitution.  The  operator  Ak  restricted  to  the  interval  I  <r  <R  leads 
to  the  Bessel  differential  equation  with  Dirichlet  boundary  conditions,  whose  solutions  are  linear  combinations 
of  Jn  and  Y„ .  This  self  adjoint  Sturm-Liouville  eigenvalue  problem  defines  a  complete  family  of  eigenfunctions 
Znm  which  are  orthogonal  with  respect  to  the  weight  function  r: 

-  ,  f  fy  /  \  /  \^  ^mp  f  dZyimir^ 

(5.3)  rZnm{i’)Znp{r)dr  —  ^ 


For  m  =  p^n  and  large  R,  the  asymptotic  expansion 


(5.4) 


Znm{r)  ~ 


sm{m{r  —  1)) 


applies  and  the  integral  (5.3)  approaches  (/?-  l)/(7rm).  Moreover,  all  eigenvalues  -m^  are  known  to  be  real 
and  strictly  negative. 

Figure  5.1  shows  typical  mode  shapes  of  $  and  the  corresponding  w  modes.  Vorticity  modes  for  low  |A:| 
exhibit  a  balance  between  vorticity  creation  at  the  wall  and  dissipation  in  the  exterior  domain.  At  high  |A:|, 
the  time  scales  m  -C  |fc|  show  no  significant  activity  in  vorticity  at  the  wall  where 


(5.5) 


-2/7r  -2(m/2)" 

Jn{m)^  +Yn{mf  r(n) 


Due  to  the  r(n)  term  in  the  denominator,  these  vorticity  modes  are  virtually  unobservable  at  the  wall.  This 
has  important  implications  for  sensing. 


Fig.  5.2.  Comparison  of  numerical  and  analytic  eigenvalues  forR  —  100  and  fc  =  0.  Logarithm  (base  10)  of  the  relative 
error  is  plotted.  Virtually  full  machine  accuracy  is  obtained  up  to  the  point  of  eigenfunction  resolution  failure  near  the  92nd 
eigenvalue.  Chebyshev  discretization  of  257  points  in  thep  direction  was  used. 
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Fig.  5.3.  Logarithm  (base  10)  of  the  numerical  eigenvalues  computed  at  various  resolutions  (thin  lines,  N  = 
17, 33, 65, 129, 257j  and  the  analytic  eigenvalues  (thick  dots).  Increasing  the  number  of  points  helps  delay  resolution  failure, 
which  still  limits  the  number  of  correctly  computed  eigenvalues. 


Comparing  the  analytic  modes  to  the  ones  computed  numerically  using  Fourier- Chebyshev  spectral 
methods,  we  achieved  nearly  full  machine  accuracy  for  properly  resolved  eigenfunctions  (figures  5.2,  5.3  and 
5.4).  As  soon  as  the  Nyquist  criterion  is  violated  anywhere  in  the  domain,  numerical  eigenvalues  exhibit 
nonphysical  behavior.  Those  modes  should  be  discarded  as  numerical  artifacts.  Care  must  be  taken  to 
construct  discretizations  which  simultaneously  resolve  all  eigenfunctions  of  interest. 

The  case  R  oo  leads  to  a  singular  Sturm-Liouville  eigenvalue  problem.  The  system  (4.23)  may  be 
diagonalized  by  means  of  Weber’s  transform  [4].  This  transform  is  a  generalization  of  the  Hankel  transform. 
The  relevant  transform  pair  applied  to  a  function  /(r),  given  our  normalization  of  Znm^  reads  as  follows: 


(5.6)  Fra  =  /(r-)Z„„(r)r  dr  W„/ 

(5.7)  f{r)=  r  Frr,Zr,,r^{r)mdm  =  W-^F  . 

Jo 

One  may  also  show  that  Weber’s  transform  satisfies  a  relation  of  the  Parseval  type: 

pOO  /*CX)  poo 

(5.8)  /  FmGrnmdm  =  /  Frnmdm  /  g{r)Znmir)r  dr 

Jo  Jo  Ji 

/oo  poo 

g{r)rdr  /  FmZnm{r)mdm 

/OO 

f{r)g{r)rdr  . 

A  simple  calculation  shows  that  this  transform  diagonalizes  Ak  so  that 


(5.11) 


VJr.AkW-^Fm^-VP?Fr 


and  the  dynamics  are  described  by 

(5.12) 

(5.13) 


Hfem  =  Wn^r) 


dt 


'km  2 
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Fig.  5.4.  Plot  density  shows  magnitude  of  numerically  obtained  modes  from  the  least  stable  (top)  to  the  most  stable 
(bottom)  for  the  circumferential  wavenumber  k  ~  I,  using  Chebyshev  discretization  in  the  p  direction.  There  are  257  grid 
points  from  the  cylinder  wall  atr  =  1  (left)  to  the  limit  of  the  computational  domain  atr  =  100  (right).  The  top  portion  is 
physically  correct.  Lower  modes  are  numeHcal  artifacts  traced  to  the  resolution  failure  starting  near  the  gridpoint  200  of  the 
eigenfunction  shown  at  vertical  coordinate  165,  which  corresponds  to  the  92nd  eigenvalue. 


One  also  obtains  a  simple  expression  for  the  pseudoenergy  of  the  flow  represented  by 

,  r  /‘27r  /.OO  oo 

(5.14)  ar,e)\dedr  =  27rl  5]|a(r)prdr 

^  — OO 

OO 

=  27r  ^  /  \Skm\^rndm  . 

k=-ooJo 

6.  Invariant  subspaces  and  optimal  control.  We  have  shown  that  the  system  dynamics  (4.23)  may 
be  diagonalized  and  that  the  invariant  subspaces  are  of  the  form  where  n  =  ||fc|  -  1|  and  m  >  0. 

Moreover,  the  pseudoenergy  form  (5.14)  is  also  simplified.  This  diagonalization  leads  to  a  family  of  optimal 
control  problems  which  preserve  the  decomposition  of  the  state  space  into  invariant  subspaces. 

We  shall  seek  to  minimize  the  time  integral  of  ||e^|p  +  ||u|p,  the  weighted  sum  of  the  pseudoenergy 
norms  of  the  flow  state  ^  and  feedback  u.  By  the  diagonalization  procediure  described  above,  we  are  led  to 
the  1-dimensional  LQR  control  problem 

(6.1)  + 
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(6.2) 


pOO 

min 77  =  m  /  \eEkm{'t)\^  +  |C^A:m(OP 

^  Jo 


dt 


where  e  >  0  is  a  small  gain  parameter,  m  >  0  is  real,  and  Ekmi't)  and  Ukmi't)  are  complex  scalars.  The 
solution  will  be  of  the  form  Ukm  =  -K{m)  Ekm-  Since  m  is  real,  the  performance  measure  reads 

_  m(e^  +  K{m)^) 


(6.3) 


where  Sfcm(O)  is  the  initial  state.  Solving  for  the  optimal  K{m)  we  obtain 
(6.4)  K{m)  =  — +  \/m^  +  e^. 

This  optimal  functional  gain  achieves  its  maximum  at  m  =  0  and  behaves  as  follows: 

e  +  O(m^)  ioT  m<^  y/e 


(6.5) 


K{m) 


6(x/2-1) 


ioY  m  =  y/e 


[  e^f{2rr?)  +  0(e'^/m®)  for  m  >  . 


The  optimal  feedback  law  in  B  representation  is  Ukmit)  =  -K{m)Ekm{t)*  The  optimal  performance  measure 
is  ry  =  mi^(m)|Hfcm(0)|^. 

6.1.  Rational  approximation  of  optimal  feedback.  While  the  above  explicit  solutions  give  the 
optimal  feedback,  their  analytic  form  is  complicated.  Fortunately,  the  performance  criterion  is  not  very 
sensitive  to  perturbations  in  the  optimal  control.  We  may  replace  the  optimal  gain  K{m)  by  a  simple 
rational  approximation  while  increasing  the  corresponding  performance  measure  rj  only  slightly. 

Consider  i^i(m)  =  c^f{2riiP‘  +  e)  as  the  rational  approximation  to  the  optimal  feedback.  This  expression 
matches  the  behavior  of  K{m)  as  m  00  and  reaches  the  same  value  at  m  =  0.  The  maximum  absolute 
error  in  gain  is  0.121585e,  reached  at  m  =  0.617541v^.  However,  the  maximum  relative  error  in  optimality 


is 


(6.6) 


7?i(m)  -7y(m) 
ry(m) 


y/e^  +  (e^  +  2  e  -f  2  m^) 


(e  H“  2  m^)  (e^  -h  e  H-  2  m"^)  (e^  +  +  m^) 


-1 


and  its  maximum  of  0.0114695  (less  than  1.147  percent)  is  reached  at  m  —  0.658911\/e.  At  other  values  of 
m  the  relative  error  in  optimality  is  even  smaller,  and  typical  relative  loss  in  performance  will  be  only  one 
percent  or  less.  By  contrast,  the  relative  performance  loss  without  any  feedback  tends  to  infinity  as  m  — >  0. 

The  advantage  of  using  the  slightly  suboptimal  gain  Ki  is  that  this  rational  function  may  be  immediately 
linked  to  the  resolvent  of  the  original  operator  (with  the  same  boundary  conditions).  The  sequence  of 
transformations  may  then  be  inverted  to  conclude  that  the  following  theorem  holds: 

Theorem  6.1.  Given  a  diagonalizable  operator  A  such  that  its  spectrum  is  real  and  negative^  and  a 
simultaneously  diagonalized  performance  measure  r]  of  type  (6.2),  the  suboptimal  feedback  kernel 


(6.7) 


ICi  =  (pel  -  2A) 


-1 


where  P  =  1  results  in  worst  case  performance  loss  of  1.14695  percent  relative  to  the  optimal  feedback 
kernel  Denoting  by  —w?  the  eigenvalues  of  A,  the  worst  case  occurs  at  m  =  0.658911v^.  The  worst  case 
relative  performance  loss  is  reduced  to  0.6437126  percent  when  p  —  0.89279  but  it  occurs  at  both  m  —  Q  and 
m  =  0.748032v^. 
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The  generalization  /S  ^  1  is  obvious,  and  its  optimal  value  reported  above  was  obtained  numerically. 
Unless  otherwise  stated,  we  shall  be  using  the  value  (3=1. 

A  better  approximation  would  express  the  optimal  gain  as  a  sum  of  different  resolvents.  The  following 
theorem  shows  how  this  can  be  done. 

Theorem  6.2.  When  z  =  and  p  is  real  and  positive^  the  expression 

approximates  yjp‘^  + 1  —  p  with  worst  case  absolute  error  of  —0.0182787  at  p  =  0.571448  and  worst  case 
relative  error  of  —3.48583  percent  at  p  =  0.854638. 


Fig.  6.1.  Exact  gain,  its  rational  approximation  (6.8)  and  their  difference  plotted  as  a  function  of  ^/p.  The  approximation 
underestimates  the  exact  value  only  slightly,  by  less  thanS.5  percent. 


It  is  useful  to  note  that  this  approximation  can  be  written  as  ^[z/{z3-p)Y  Further  calculations  demon¬ 
strate  that  the  worst  case  relative  performance  loss  is  0.0256434  percent,  reached  at  p  =  0.617808  (figure  6.1). 
Taking  p  to  be  the  eigenvalue  of  —Aje.,  we  obtain  the  corresponding  approximate  feedback  kernel  as  a  sum 
of  two  Green’s  functions: 


(6.9) 

(6.10) 


]C2  =  ^{ez{ezT  -  A)  ^ ez{ezl  -  A) 


=  e3? 


ez  {ezX  —  A) 


-1 


These  results  will  be  used  to  derive  analytic  approximations  to  the  optimal  feedback  operator. 


6.1.1.  A  simple  example  of  rational  approximation.  The  heat  equation  on  an  infinite  line  is  ht  = 
hxx  and  can  be  diagonalized  by  Fourier  transform  so  that  ht  =  —k^h.  For  this  simple  example,  the  spatial 
representation  of  the  optimal  feedback  kernel  can  be  obtained  analytically  [16]  in  terms  of  the  generalized 
hypergeometric  function  qFs.  The  exact  expression  for  the  optimal  control  is  u{x)  =  ~/ K{x,x^)h{x')dx^ 
where 


k{x,x^)  =  -  (c/tt)^^^ 


1  3 

2’4’4’ 
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(6.11)  |s|7r^/^  0F3  j,2; 

where  we  have  defined  s  =  {x  -  x')^/ef4.  The  two  rational  approximations  defined  above  evaluate  to 

(6.12)  Ki(x,x')  =  (e/2)^/^e“^'^l®l 
and 

I 

-|s|2> 


(6.13) 


K2{x,  x')  =  J ^  e  1*1  cos  -  |s|  2^JV2-l^  . 


Fig.  6.2.  Spatial  form  of  the  optimal  feedback  kernel  for  the  heat  equation  on  an  infinite  line.  Left:  the  exact  k  (6.11) 
and  its  rational  approximations  Ki  and  K2  as  functions  of  s  for  e  =  1.  Right:  error  in  the  two  rational  approximations  shown 
using  expanded  vertical  scale.  See  text. 


For  both  rational  approximations  the  maximum  error  occurs  at  s  =  0  (figure  6.2).  The  exact  value 
of  the  feedback  kernel  at  s  ==  0  should  be  -(e/7r)^/^r(-3/4)r(5/4)/2.  By  comparison,  ki  gives  a  10.14 
percent  lower  value  while  k>2  underestimates  by  only  1.27  percent.  We  have  already  shown  that  both  of 
these  suboptimal  strategies  perform  nearly  as  well  as  the  optimal  control.  In  spectral  representation,  the 
corresponding  worst  case  relative  performance  losses  are  less  than  1.147  percent  and  less  than  0.026  percent, 
respectively. 

We  note  that  the  k,2  exhibits  sign  changes  starting  at  s  =  57r\/\/2H-  1/16  «  1.52541.  This  is  qualitatively 
correct  since  the  exact  feedback  kernel  k  also  exhibits  sign  changes  (spaced  by  approximately  1.12  at  first) 
starting  at  s  »  1.47894.  The  magnitude  of  the  feedback  kernel  beyond  that  point  is  at  most  0.000334e^/^  so 
this  interesting  behavior  has  little  bearing  on  the  form  of  the  optimal  feedback. 


6.1.2,  Rational  approximation  and  2D  diffusion.  Consider  the  heat  equation  in  an  infinite  plane: 
ht  =  Ah.  While  the  exact  spatial  form  of  the  optimal  feedback  kernel  could  not  be  obtained  analytically, 
rotational  symmetry  leads  to  the  rational  approximations 


(6.14) 


«;i(x,x')  =  —Ko 


(’■4) 


and 

(6.15) 


l  +  i 


Ko 


/(I  +  i)e 


l-i 

H — tt-'Ko 


where  the  distance  between  the  points  in  the  plane  x  and  x'  defines  r  =  |x  —  x'|.  These  expressions  involve 
the  modified  Bessel  function  Kq^  which  has  a  logarithmic  singularity  at  r  =  0.  The  difference  between  these 
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two  approximations  approaches  its  maximum  as  r  tends  to  zero: 


(6.16) 


lim  («;2  —  m)  = 


e^(7r  -  21og(2)) 
327r 

0.01746036^  . 


6.1.3.  Rational  approximation  and  3D  diffusion.  For  the  heat  equation  in  3D  free  space,  the 
rational  approximations  yield 


(6.17) 
and 

(6.18) 


^  -Ty/e/2 


«;2(x,x') 


r-v/^(I+^/2 


COS 


1) 


4\/27rr 

where  r  =  |x  —  x'|  is  the  distance  between  the  two  points  in  space.  The  origin  r  =  0  is  now  a  simple  pole, 
where  the  difference  between  these  two  approximations  approaches  its  maximum  given  by 

6^/2  +  7^2 -ij 


(6.19) 


lim  {k2  —  /^i)  = 

r— >0 


8\/27r 

0.04624236^/2 


7.  Approximately  optimal  control  of  Stokes  flow.  To  obtain  the  rational  approximation  of  the 
optimal  feedback  kernel  for  equation  (4.23),  we  first  need  to  determine  the  Green’s  function  which  solves  the 
auxiliary  problem 

(7.1)  ^nf^k  “  ^^^k 

where  n  =  ||fc|  —  1|  and  Dirichlet  boundary  conditions  are  applied. 

The  required  Green’s  function  can  be  derived  explicitly  in  terms  of  modified  Bessel  functions 

(7.2)  Gn{r,s)  -  (/„(ainin(r,s))/!r„(amax(r,s))  -  ^ 

where  a  =  yJTz,  The  solution  can  now  be  written  as 

/oo 

Gn{r,s)i]^{s)sds 

where  /x/c  and  are  now  expressed  as  functions  of  the  radius.  The  rational  approximations  to  the  optimal 
feedback  kernel  are  obtained  in  terms  of  Gn- 

From  now  on,  we  shall  focus  on  the  case  z  =  y/±if2  and  refer  to  the  resulting  approximate  kernel  of 
type  (6.9)  as  simply  k. 

Theorem  7.1.  The  rational  approximation  of  type  (6.9)  to  the  optimal  feedback  kernel  for  the  system 
governed  by  H.2S)  is  given  by 

(7.4)  kk{r,  s)  =  eSR  [Gn  (r,  s)]  using  a  = 

where  n  =  ||fc|  —  1|  and  Gn(r,  s)  is  defined  by  (7.2)  for  the  indicated  value  of  a. 

These  feedback  kernels  exhibit  a  pronounced  ridge  along  the  line  r  =  s  representing  colocated  sensing  and 
actuation.  Figure  7.1  shows  the  typical  form  of  control  laws  derived  from  feedback  kernels  kk-  The  colocated 
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Fig.  7.1.  Feedback  at  r  due  to  disturbance  s  for  e  =  0.001  and  fc  =  1, . . . ,  4. 


control  ridge  becomes  essentially  flat  for  large  |A:|.  However,  this  feedback  gain  does  not  yet  provide  spatial 
clues  because  it  is  expressed  in  our  nonlocal  pseudodifferential  representation.  The  spatial  information  we 
seek  involves  returning  to  the  vorticity  representation.  This  can  be  done  as  follows. 

We  shall  assume  that  the  gain  parameter  e  is  small  and  consider  disturbances  whose  support  is 
contained  within  a  finite  computational  domain  of  size  R  <  l/|a|.  Under  these  conditions,  |amax(l,r,  s)| 

1  and  the  Green’s  function  behaves  as 


(7.5) 


Gn{r,s) 


G?  (min(r,  s)”"  -  min(r,  s)  ”) 


2n  max(r, 

when  n  7^  0.  This  applies  when  /  ±1  and  we  obtain 


(7.6) 


/^fc(r,  s) 


^  (min(r,  s)”  —  min(r,  s)  ^)  def 


4n  max(r,  s)’^ 


=  A,(r,s). 


For  the  special  case  n  =  0,  we  obtain  the  expression 


(7.7) 


Go(r,  s)  ~  c?  log(min(r,  s)) 


/  log(max(r,s))\ 
V  7  +  log(a/2)  ) 


where  7  0.577216  is  Euler’s  constant  gamma.  The  term  log(a/2)  evaluates  to  log(e)/2  ~  51og(2)/4±i7r/8 

when  z  =  y/±i/2.  After  some  algebra,  we  obtain 


(7.8) 


6^  log(min(r,  s)) 
2 


(1  +  a  log(max(r,  s)))  A±i(r,  s) 


where  a  is  defined  by  (7.13). 

The  asymptotic  form  of  h  will  constitute  our  low  gain  rational  approximation  (denoted  by  C  and  con¬ 
sidered  an  integral  operator  with  kernel  A(r,  s))  which  is  applicable  within  the  domain  of  radius  R  <C  l/\/e- 
The  approximately  optimal  control  in  pseudodifferential  representation  is  — but  to  get  its  spatial  inter¬ 
pretation  we  must  go  back  to  the  vorticity  representation,  where  Ck  is  transformed  to 


(7.9) 


Ck  =  'R  , 
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Theorem  7.2.  The  low  gain  rational  approximation  to  the  optimal  feedback  kernel  in  vorticity  repre¬ 
sentation  is  given  by 


Ao(r,  s)  =  log  (miii(r,  s)/s)  /2 
AiiCr,  s)  =  (2 inin(r,  s)^  —  s^{2a log(s)  -  a  +  2)  —  a)  / (8rs) 

>^kir,s)  =  ^inin(r,s)^l''l  -  |fc|s^  +  |fc|  “  l)  /  (4|fc|(rs)l''l  j 

>  2  and  the  parameter  a  is  given  by 

^  8(7r  +  87  +  21og(£V32)) 

^  TT^  +  (87  +  21og(e2/32))2' 

This  approximation  applies  within  the  domain  of  radius  R  <C 

We  note  that  if  the  support  of  ujk  is  contained  within  a  computational  domain  of  radius  i?,  so  is  the 
support  of  and  therefore  the  low  gain  approximation  remains  valid.  One  must  also  make  a  distinction 
between  |A:|  —  1  and  n  =  ||fc|  —  1|,  and  pay  attention  to  the  limits  of  integration.  Once  that  is  done,  the 
theorem  follows  by  a  straightforward  but  lengthy  calculation. 

This  result  is  given  as  a  function  of  radius  and  the  nearly  optimal  control  law  now  reads 

/oo 

\k{r,s)u)k{s)sds  . 

As  we  have  shown  before,  the  maximum  relative  error  of  the  rational  approximation  /C2  is  under  3.5  percent, 
resulting  in  maximum  relative  performance  loss  of  under  0.026  percent.  The  kernel  \k{r,  s)  also  involves  the 
low  gain  approximation,  where  the  neglected  term  is  of  relative  order  0(e).  Since  we  expect  that  e  is  in  the 
range  of  1  percent  or  less,  Afc(r,  s)  is  an  excellent  approximation  to  the  exact  optimal  distributed  feedback. 
Using  e  =  0.001  and  R  =  10,  the  maximum  relative  error  of  the  low  gain  approximation  remains  under  2 
percent. 


(7.10) 

(7.11) 

(7.12) 

where  \k\ 

(7.13) 


7.1.  How  to  control  a  particulair  vortex.  We  shall  assume  that  a  resolution  limit  on  control  and 
observation  exists,  so  that  only  the  subsystems  \k\  <  K  need  to  be  considered.  This  assumption  suggests  an 
evaluation  of  the  control  effort  at  each  wavenumber. 

Given  a  vorticity  distribution  uj  which  satisfies  integral  compatibility  constraints,  we  shall  focus  at 
the  contribution  of  the  vorticity  near  {r,6)  =  {s.cj)).  Consider  a  particular  vortex  of  the  form  uj{r,0)  — 
d{r  —  s)6{r{9  —  (f>))  whose  Fourier  transform  yields 


(7.15) 


For  simplicity,  we  shall  take  0  =  0  and  consider  the  distributed  feedback  produced  by  the  low  gain 
rational  approximation.  This  single  vortex  in  the  flow  induces  a  vortex  sheet  at  the  solid  boundary,  which 
does  not  contribute  to  the  control  because  the  feedback  kernel  is  zero  for  vortices  at  the  boundary. 

Distributed  vorticity  control  (Fig.  7.2)  exhibits  two  pronounced  ridges  within  the  bounded  domain  where 
our  low  gain  rational  approximation  applies. 

The  ridge  along  the  line  r  =  s  represents  the  colocated  control.  Within  the  region  of  radius  i?  <C  l/\/6  this 
control  is  zero  for  fc  =  0,  grows  logarithmically  with  r  for  |fc|  =  1,  and  saturates  at  a  fixed  value  proportional 
to  l/\k\  for  \k\  >  2.  This  control  represents  colocated  vorticity  damping,  opposing  the  disturbance  vortex. 

Another  ridge  is  present  along  the  solid  boundary.  The  magnitude  of  this  control  grows  logarithmically 
with  s  for  fc  =  0,  behaves  as  slog(s)  for  |fc|  =  1  and  as  s“l^l(s^  -  1)  for  |fcl  >  2.  This  form  of  control  creates 
a  vortex  at  the  wall  of  the  same  orientation  as  the  disturbance  vortex. 
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Fig.  7.2.  Distributed  vorticity  control  atr  in  response  to  a  vortex  at  s  for  e  =  0.001  and  A;  =  0, ...  ,5,  Boundary  control 
and  colocated  control  dominate. 


Fig.  7.3.  Left:  Boundary  control  (positive  traces)  and  colocated  control  (negative  traces)  as  functions  of  the  radial  location 
of  the  disturbance  vortex  for  c  =  0.001  and  k  =  0, . . . ,  5.  Right:  Sum  of  boundary  control  and  colocated  control  Within  the 
region  R  boundary  control  dominates  for\k\  <2  or  sujficiently  close  to  the  wall. 


Figure  7.3  demonstrates  that  for  |fc|  <  2  more  control  is  applied  at  the  boundary  than  at  the  vortex 
location  (within  the  region  where  low  gain  approximation  holds).  For  |fc|  >  3,  boundary  control  still  domi¬ 
nates  within  a  thin  layer  along  the  boundary,  but  colocated  control  is  preferred  in  the  rest  of  the  domain. 
Using  the  low  gain  approximation,  analysis  and  numerical  experiments  indicate  that  for  |A;|  >  3  the  region 
of  boundary  control  dominance  does  not  depend  on  e  and  extends  to 


(7.16) 


s=p  +  q^ 


1 1  —  -h  \k\{p^  —  1) 

|fci(|fc|-2) 


where  p  =  and  ^  «  0.9  ±  0.01. 

Although  boundary  control  dominates  for  low  |fc|  and  near  the  solid  boundary,  the  control  input  depends 
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on  the  vorticity  distribution  within  the  flow.  Vorticity  at  the  wall  has  zero  impact  on  the  control.  Therefore, 
sensing  vorticity  at  the  wall  is  not  directly  useful.  One  must  estimate  the  flow  state  away  from  the  wall. 

8.  Conclusions.  The  analytic  expressions  we  obtained  were  used  to  validate  a  numerical  scheme  based 
on  Fourier-Chebyshev  spectral  methods.  While  full  accuracy  was  achieved  for  properly  resolved  eigenfunc¬ 
tions,  numerical  artifacts  were  observed  in  other  cases.  This  reinforces  the  obvious  but  important  point  that 
the  chosen  discretization  must  simultaneously  resolve  all  eigenfunctions  of  interest  throughout  the  computa¬ 
tional  domain. 

The  analytic  form  of  the  approximately  optimal  feedback  shows  that  boundary  control  dominates  for 
|fc|  <  2  or  when  the  vorticity  disturbance  is  sufficiently  close  to  the  wall.  We  obtained  explicit  estimates  on 
the  distance  from  the  wail  within  which  the  optimal  control  favors  actuation  by  means  of  boundary  vortex 
generators. 

We  have  also  shown  that  vorticity  at  the  wall  does  not  influence  the  optimal  control.  The  control  effort 
depends  on  estimating  the  flow  state  away  from  the  wall.  This  is  a  dual  problem  to  the  results  reported 
here.  The  mode  shapes  in  vorticity  representation  show  that  some  modes  are  virtually  unobservable  at  the 
wall,  leading  to  limits  on  the  effectiveness  of  boundary  sensing.  In  future  work,  we  hope  to  report  on  the 
estimation  problem,  and  to  extend  this  approach  to  more  complex  flows. 
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